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Abstract. We have studied quantization and confinement effects in nanostructured 
superconductors. Three different types of nanostructured samples were investigated: 
individual structures (line, loop, dot), 1-dimcnsional (ID) clusters of loops and 2D 
clusters of antidots, and finally large lattices of antidots. Hereby, a crossover from 
individual elementary " plaquettes" , via clusters, to huge arrays of these elements, 
is realized. The main idea of our study was to vary the boundary conditions for 
confinement of the superconducting condensate by taking samples of different topology 
and, through that, modifying the lowest Landau level Elll{H). Since the critical 
temperature versus applied magnetic field T C (H) is, in fact, Elll(H) measured in 
temperature units, it is varied as well when the sample topology is changed through 
nanostructuring. We demonstrate that in all studied nanostructured superconductors 
the shape of the T C (H) phase boundary is determined by the confinement topology in 
a unique way. 



I INTRODUCTION 

A Confinement and Quantization 

"Confinement" and "quantization" are two closely interrelated definitions: if a 
particle is "confined" then its energy is "quantized", and vice versa. According to 
the dictionary, to "confine" means to "restrict within limits" , to "enclose" , and even 
to " imprison" . A typical example, illustrating the relation between confinement and 
quantization, is the restriction of the motion of a particle by an infinite potential 
well of size La- Due to the presence of an infinite potential U(x) (Fig. 1) for x < 
and x > La-, the wave function ty(x) describing the particle is zero outside the well: 
\& = for x < and x > La and, in the region with U(x) =0 (0 < x < La), 
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the solutions of the one-dimensional Schrodinger equation correspond to standing 
waves with an integer number n of half wavelengths A along La '■ n A„/2 = La- 
This simple constraint results in the well-known quantized energy spectrum 

U -> 00 JJ (JO 




FIGURE 1. Confinement and quantization of the motion of a particle by an infinite potential 
well with size La for n = 1, 2 and 3. 
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Here k n is the wave number and m is the free electron mass. To have an idea about 
the characteristic energy scales involved and their dependence upon the confinement 
length La-, we have calculated (see Table 1) the energies E 1 (Eq. (1)) for electrons 
confined by an infinite potential well with the sizes 1 A, 1 nm and 1 jim. 



TABLE 1. Confinement by the infinite potential well 
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Temperature T 
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B Nanostructuring 

Recent impressive progress in nanofabrication has made it possible to realize the 
whole range of confinement lengths La '■ from \\im (photo-and e-beam lithography), 
via 1 nm to 1 A (single atom manipulation) and, through that, to control the 



confinement energy (temperature) from a few mK higher up to far above room 
temperature (Table 1). 

This progress has stimulated dramatically the experimental and theoretical stud- 
ies of different nanostructured materials and individual nanostructures. The in- 
terest towards such structures arises from the remarkable principle of "quantum 
design", when quantum mechanics can be efficiently used to tailor the physical 
properties of nanostructured materials. 

Nanostructuring can also be considered as a sort of artificial modulation. We can 
identify then the main classes of nanostructured materials using the idea of their 
modification along one-, two- or three-axes, thus introducing 1-dimensional (ID)-, 
2D- or 3D- artificial modulation (Fig. 2). The ID or "vertical" modulation repre- 



FIGURE 2. Schematic presentation of the vertical modulation in superlattices or multilayers 
(a), of the horizontal modulation achieved by a lateral repetition of elements A and B (b) and of 
the 1D+2D=3D artificial modulation (c). 

sents then the class of superlattices or multilayers (Fig. 2a) formed by alternating 
individual films of two (A, B) or more different materials on a stack. Some examples 
of different types of multilayers are superconductor /insulator (Pb/Ge, WGe/Ge,...), 
superconductor/metal (V/Ag,...), superconductor/ferromagnet (Nb/Fe, V/Fe,...), 
ferromagnet /metal (Fe/Cr, Cu/Co,...), etc. 

The " horizontal" (lateral) superlattices (Fig. 2b) correspond to the 2D artificial 
modulation achieved by a lateral repetition of one (A), two (A,B) or more ele- 
ments. As examples, we should mention here antidot arrays or antidot lattices, 
when A=microhole ("antidot"), or arrays and lateral superlattices consisting of 
magnetic dots. 

If the 2D lateral modulation is applied to each individual layer of a multilayer 
or superlattice, then we deal with the 1D+2D=3D artificial modulation (Fig. 2c). 
For example, if arrays of antidots are made in a multilayer, then we have a system 
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with 3D artificial modulation which combines 2D lateral "horizontal" with the ID 
"vertical" modulation. 

Finally, macroscopic nanostructured samples, with a huge number N of repeated 
elementary "plaquettes" (A,B,...), are examples of very complicated systems if the 
confined charge carriers or flux lines are strongly interacting with each other and 
the relevant interaction is of a long range. In this case the essential physics of 
such systems can be understood much better if we use clusters of elements (N ~ 
10), instead of their huge arrays (N — > oo) (Fig. 3). These clusters, occupying an 



A 


A 




A 


A 


A 


A 






A 


A 


A 


A 


A 


A 




A 


A 


A 


A 


A 


A 




A 


A 


A 


A 


CLUSTER 


HUGE ARRAY 



INDIVIDUAL 
STRUCTURE 

FIGURE 3. Schematic presentation of an individual structure, a cluster and a huge array. 



intermediate place between individual nanostructures (N — 1) and nanostructured 
materials (N — > oo), are very helpful model objects to study the interactions be- 
tween flux lines or charge carriers confined by elements A. The "growth" of clusters 
on the way from an individual object A to a huge array of A's can be done either in 
a ID or 2D fashion (Fig. 3), thus realizing ID chains or 2D-like clusters of elements 
A. 



C Confining the Superconducting Condensate 

The nanostructured materials and individual nanostructures, introduced in the 
previous section, can be prepared using the modern facilities for nanofabrication. 
It is worth, however, first asking ourselves a few simple questions like: why do we 
want to make such structures, what interesting new physics do we expect, and why 
do we want to focus on superconducting (and not, for example, normal metallic) 
nanostructured materials? 

First of all, by making nanostructured materials, one creates an artificial poten- 
tial in which charge carriers or flux lines are confined. The size La of an elementary 
"plaquette" A, gives roughly the expected energy scale in accordance with Table 1, 
while the positions of the elements A determine the pattern of the potential mod- 
ulation. The concentration of charge carriers or flux lines can be controlled by 
varying a gate voltage (in 2D electron gas systems) [1] or the applied magnetic 
field (in superconductors) [2]. In this situation, different commensurability effects 
between the fixed number of elements A in an array and a tunable number of charge 
or flux carriers are observed. 



Secondly, modifying the sample topology in nanostructured materials creates a 
unique possibility to impose the desired boundary conditions, and thus almost "im- 
pose" the properties of the sample. A Fermi liquid or a superconducting condensate 
confined within such materials will be subjected to severe constraints and, as a re- 
sult, the properties of these materials will be strongly affected by the boundary 
conditions. 

While a normal metallic system should be considered quantum-mechanically by 
solving the Schrodinger equation: 

— {-ihV -eA) 2 m + im = , 
2m v ' 

a superconducting system is described by the two coupled Ginzburg-Landau (GL) 
equations: 

^_(_^V-eM) 2 ^ s + /3|^ s | 2 ^ = -a^ (3) 



(2) 



j s = V x h = — [V* a (-ihV - eM)f s + * s (jW - e*A)V* a \ , (4) 
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with A the vector potential which corresponds to the microscopic field h = rot A/ fj, , 
U the potential energy, E the total energy, a a temperature dependent parameter 
changing sign from a > to a < as T is decreased through T c , f3 a positive 
temperature independent constant, m* the effective mass which can be chosen 
arbitrarily and is generally taken as twice the free electron mass m. 

Note that the first GL equation (Eq. (3)), with the nonlinear term /3|\l/ s | 2 \l/ s 
neglected, is the analogue of the Schrodinger equation (Eq. (2)) with U — 0, when 
making a few substitutions: <-> ^, e* <-> e, —a <-> E and m* <-> m. The 
superconducting order parameter ^ s corresponds to the wave function \1/ in Eq. (2). 
The effective charge e* in the GL equations is 2e, i.e. the charge of a Cooper pair, 
while the temperature dependent GL parameter a 

~ a = 2^aT) (5) 

plays the role of E in Schrodinger equation. Here £(T) is the temperature dependent 
coherence length: 

t(T) = -4^=- (6) 
V 1 ~ tTo 

The boundary 
conditions for interfaces normal metal-vacuum and superconductor-vacuum are, 
however, different (Fig. 4): 
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FIGURE 4. Boundary conditions for interfaces normal metal-vacuum and 
superconductor- vacuum. 



i.e. for normal metallic systems the density is zero, while for superconducting sys- 
tems, the gradient of \I/ S (for the case A = 0) has no component perpendicular to the 
boundary. As a consequence, the supercurrent cannot flow through the boundary. 
The nucleation of the superconducting condensate is favored at the superconduc- 
tor/ vacuum interfaces, thus leading to the appearance of superconductivity in a 
surface sheet with a thickness £(T) at the third critical field H c3 (T). 

For bulk superconductors the surface-to-volume ratio is negligible and there- 
fore superconductivity in the bulk is not affected by a thin superconducting sur- 
face layer. For nanostructured superconductors with antidot arrays, however, 
the boundary conditions (Eq. (8)) and the surface superconductivity introduced 
through them, become very important if La < £,(T). The advantage of supercon- 
ducting materials in this case is that it is not even necessary to go to nm scale (like 
for normal metals), since for La of the order of 0.1-1.0 fim the temperature range 
where La < £(T), spreads over 0.01 — 0.1 K below T c due to the divergence of £(T) 



In principle, the mesoscopic regime La < £,{T) can be reached even in bulk 
superconducting samples with La ~ 1 cm — 1 m, since £(T) diverges. However, the 
temperature window where La < £(T) is so narrow, not more than ~ 1 nK below 
T c0 , that one needs ideal sample homogeneity and perfect temperature stability. 

In the mesoscopic regime L A < £CO, which is quite easily realized in (perforated) 
nanostructured materials, the surface superconductivity can cover the whole avail- 
able space occupied by the material, thus spreading superconductivity all over the 
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at T -> T c0 (Eq. (6)). 



sample. It is then evident that in this case surface effects play the role of bulk 
effects. 

Using the similarity between the linearized GL equation (Eq. (3)) and the 
Schrodinger equation (Eq. (2)), we can formalize our approach as follows: since 
the parameter -a (Eqs. (3) and (5)) plays the role of energy E (Eq. (2)), then the 
highest possible temperature T C (H) for the nucleation of the superconducting state 
in presence of the magnetic field H always corresponds to the lowest Landau level 
Elll(H) found by solving the Schrodinger equation (Eq. (2)) with "superconduct- 
ing" boundary conditions (Eq. (8)). 
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FIGURE 5. Landau level scheme for a particle in a magnetic field. From the lowest Landau 
level Elll(H) the second critical field H C 2{T) is derived (solid line). 

Figure 5 illustrates the application of this rule to the calculation of the upper 
critical field H c2 (T): indeed, if we take the well-known classical Landau solution 
for the lowest level in bulk samples E LLL (H) = huj/2, where u = e*fXoH/m* is the 
cyclotron frequency. Then, from -a = Elll(H) we have 



ti 2 



2m* £ 2 (T) 

and with the help of Eq. (5), we obtain 

H H c2 (T) = 



huj 
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with $ = h/e* = h/2e the superconducting flux quantum. 

In nanostructured superconductors, where the boundary conditions (Eq. (8)) 
strongly influence the Landau level scheme, E LLL (H) has to be calculated for each 



different confinement geometry. By measuring the shift of the critical temperature 
T C (H) in a magnetic field, we can compare the experimental T C (H) with the cal- 
culated level Elll(H) and thus check the effect of the confinement topology on 
the superconducting phase boundary for a series of nanostructured superconduct- 
ing samples. The transition between normal and superconducting states is usually 
very sharp and therefore the lowest Landau level can be easily traced as a function 
of applied magnetic field. Except when stated explicitly we have taken the mid- 
point of the resistive transition from the superconducting to the normal state, as 
the criterion to determine T C (H). 

In this paper we present the systematic study of the influence of the confinement 
geometry on the superconducting phase boundary T C {H) in a series of nanostruc- 
tured samples. We begin with individual nanostructures of different topologies 
(lines, loops, dots) (Section II) and then focus on "intermediate" systems: clusters 
of loops fabricated in the form of a ID chain of loops (Section III. A) or 2D antidot 
clusters (Section III.B). Finally, we move on to huge arrays of antidots in Section 
IV where we report on the T C (H) boundary for superconducting films with antidot 
lattices. The main emphasis of the paper is on the demonstration of the importance 
of the confinement geometry for the superconducting condensate and on the related 
quantization phenomena in nanostructured superconductors through the study of 
the phase boundaries T C (H). 



II INDIVIDUAL STRUCTURES: LINE, LOOP AND 



We begin this section by presenting the experimental results on the T C (H) phase 
boundary of individual superconducting mesoscopic structures of different topology. 
Simultaneously, we have kept other parameters of these samples, like material from 
which they are made (Al), the width of the lines (w = 0.15 fim) and the film 
thickness t = 25 nm the same for all three structures, thus directly relating the 
differences in T C (H) to topological effects. The magnetic field H is always applied 
perpendicular to the structures. 



In Fig. 6a the phase boundary T C (H) of a mesoscopic line is shown. The solid 
line gives the T C (H) calculated from the well-known formula [3]: 
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which, in fact, describes the parabolic shape of T C (H) for a thin film of thickness w 
in parallel magnetic field. Since the cross- sect ion, exposed to the applied magnetic 



field, is the same for a film of thickness to in a parallel magnetic field and for a 
mesoscopic line of width w in a perpendicular field, the same formula can be used 
for both [4]. Indeed, the solid line in Fig 6a is a parabolic fit of the experimental 
data with Eq. (11) where £(0) = 110 nm was obtained as a fitting parameter. The 
coherence length obtained using this method, coincides reasonably well with the 
dirty limit value f (0) = 0.85(£(/) 1/2 = 132 nm calculated from the known BCS 
coherence length £ — 1600 nm for bulk Al [5] and the mean free path £ — 15 nm, 
estimated from the normal state resistivity p at 4.2 K [6]. 

We can use also another simple argument to explain the parabolic relation 
T C (H) oc H 2 : the expansion of the energy E(H) in powers of H, as given by 
the perturbation theory, is [7]: 

E{H) = E + A X LH + A 2 SH 2 + ■■■ (12) 

where A\ and A 2 are constant coefficients, the first term Eq represents the energy 
levels in zero field, the second term is the linear field splitting with the orbital 
quantum number L and the third term is the diamagnetic shift with S, being the 
area exposed to the applied field. 

Now, for the topology of the line with a width w much smaller than the Larmor 
radius r H ^> w, any orbital motion is impossible due to the constraints imposed by 
the boundaries onto the electrons inside the line. Therefore, in this particular case 
L = and E(H) = EQ + A2SH 2 , which immediately leads to the parabolic relation 
T c oc H 2 . This diamagnetic shift of T C (H) can be understood in terms of a partial 
screening of the magnetic field H due to the non-zero width of the line [8]. 




FIGURE 6. The measured superconducting/normal phase boundary as a function of the re- 
duced temperature T C (H)/T c0 for a) the line, and b) the loop and the dot. The solid line in (a) is 
calculated using Eq. (11) with £(0) = 110 nm as a fitting parameter. The dashed line represents 
T C (H) for bulk Al. 



B Loop 



The T C (H) of the mesoscopic loop, shown in Fig. 6b, demonstrates very distinct 
Little-Parks (LP) oscillations [9] superimposed on a monotonic background. A 
closer investigation leads to the conclusion that this background is very well de- 
scribed by the same parabolic dependence as the one which we just discussed for the 
mesoscopic line [4] (see the solid line in Fig. 6a). As long as the width of the strips 
w, forming the loop, is much smaller than the loop size, the total shift of T C (H) 
can be written as the sum of an oscillatory part, and the monotonic background 
given by Eq. (11) [4,10]: 



T C {H) = T c 
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where R 2 = Ri R 2 is the product of inner and outer loop radius, and the magnetic 
flux threading the loop <3> = t^R 2 ^qH. The integer n has to be chosen so as to 
maximize T C (H) or, in other words, selecting Elll(H). 

The LP oscillations originate from the ffuxoid quantization requirement, which 
states that the complex order parameter \l/ s should be a single- valued function when 
integrating along a closed contour 

Vy? ■dl = n1-n n = ■ ■ ■ , -2, -1, 0, 1, 2, • • • (14) 

where we have introduced the order parameter as = |\I/ s |exp (up). Fluxoid 
quantization gives rise to a circulating supercurrent in the loop when $ ^ n&o, 
which is periodic with the applied flux $/$o- 

Using the sample dimensions and the value for £(0) obtained before for the meso- 
scopic line (with the same width w = 0.15 /im), the T C (H) for the loop can be 
calculated from Eq. (13) without any free parameter. The solid line in Fig. 6b 
shows indeed a very good agreement with the experimental data. It is worth not- 
ing here that the amplitude of the LP oscillations is about a few mK - in qualitative 
agreement with the simple estimate given in Table 1 for La — 1 fim. 

Another interesting feature of a mesoscopic loop and other mesoscopic structures 
is the unique possibility they offer for studying nonlocal effects [11]. In fact, a single 
loop can be considered as a 2D artificial quantum orbit with a fixed radius, in 
contrast to Bohr's description of atomic orbitals. In the latter case the stable radii 
are found from the quasiclassical quantization rule, stating that only an integer 
number of wavelengths can be set along the circumference of the allowed orbits. 
For a superconducting loop, however, supercurrents must flow, in order to fulfil the 
fluxoid quantization requirement (Eq. (14)), thus causing oscillations of the critical 
temperature T c versus magnetic field H. 

In order to measure the resistance of a mesoscopic loop, electrical contacts have, 
of course, to be attached to it, and as a consequence the confinement geometry is 
changed. A loop with attached contacts and the same loop without any contacts 



are, strictly speaking, different mesoscopic systems. This "disturbing" or "invasive" 
aspect can now be exploited for the study of nonlocal effects [11]. Due to the 
divergence of the coherence length £(T) at T = T c0 (Eq. (6)) the coupling of the 
loop with the attached leads is expected to be very strong for T — > T c0 . 
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FIGURE 7. Local (Vi/V 2 ) and nonlocal phase boundaries (V1/V3 or V 2 /V 4 ) T C (H). The mcsur- 
ing current is sent through h/h- The solid and dashed lines correspond to the theoretical T C {H) 
of an isolated loop and a one-dimensional line, respectively. The inset shows a schematic of the 
structure, where the distance P = 0.4 \im. 

Fig. 7 shows the results of these measurements. Both "local" (potential probes 
across the loop Vi/V 2 ) and "nonlocal" (potential probes aside of the loop V1/V3 
or V2/V4) LP oscillations are clearly observed. For the "local" probes there is an 
unexpected and pronounced increase of the oscillation amplitude with increasing 
field, in disagreement with previous measurements on Al microcylinders [10]. In 
contrast to this, for the "nonlocal" LP effect, the oscillations rapidly vanish when 
the magnetic field is increased. 

When increasing the field, the background suppression of T c (Eq. (11)) results in 
a decrease of £(T). Hence, the change of the oscillation amplitude with H is directly 
related to the temperature-dependent coherence length. As long as the coherence 
of the superconducting condensate protuberates over the nonlocal voltage probes, 
the nonlocal LP oscillations can be observed. 

On the other hand, the importance of an "arm" attached to a mesoscopic loop, 
was already demonstrated theoretically by de Gennes in 1981 [12]. For a perfect ID 
loop (vanishing width of the strips) adding an "arm" will result in a decrease of the 
LP oscillation amplitude, what we observed indeed at low magnetic fields, where 
£(T) is still large. With these experiments, we have proved that adding probes 




to a structure considerably changes both the confinement topology and the phase 
boundary T C (H). 



C Dot 

The Landau level scheme for a cylindrical dot with "superconducting" boundary 
conditions (Eq. (8)) is presented in Fig. 8. Each level is characterized by a certain 
orbital quantum number L where \1/ s = |\I/ s |exp (^%L(p) [13]. The levels, corre- 
sponding to the sign " +" in the argument of the exponent are not shown since they 
are situated at energies higher than the ones with the sign The lowest Landau 
level in Fig. 8 represents a cusp-like envelope, switching between different L values 
with changing magnetic field. Following our main guideline that Elll(H) deter- 
mines T C (H), we expect for the dot the cusp-like superconducting phase boundary 
with nearly perfect linear background. The measured phase boundary T C (H) , shown 
in Fig. 6b, can be nicely fitted by the calculated one (Fig. 8), thus proving that 
T C (H) of a superconducting dot indeed consists of cusps with different L's [14]. 
Each fixed L describes a giant vortex state which carries L flux quanta $o- The 
linear background of the T C (H) dependence is very close to the third critical field 
Hcs(T) ~ 1.69 H c2 (T) [15]. Contrary to the loop, where the LP oscillations are 
perfectly periodic, the dot demonstrates a certain aperiodicity [16], in very good 
agreement with the theoretical calculations [14,17]. 




FIGURE 8. Energy level scheme versus normalized flux $/<f>o for a superconducting cylinder 
in a parallel magnetic field. The cusp-like H c3 (T) line is formed due to the change of the orbital 
quantum number L. 



Ill CLUSTERS OF LOOPS AND ANTIDOTS 
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FIGURE 9. AFM micrographs of the studied structures: a) the bola, b) the double loop, and 
c) the triple loop. 

After the description of the confinement effects of several individual supercon- 
ducting structures (A = line, loop, dot) we are ready to move further on to clusters 
of elements A (Fig. 3) on our way from "single plaquette" samples to materials 
nanostructured by introducing huge arrays of elements A. First we take A = loop 
and consider one-dimensional multiloop structures: "bola", double and triple loop 
Al structures. Figure 9 shows a AFM image of the structures. For these geometries 
some interesting theoretical predictions have been made, for which no experimental 
verification has been carried out up to know (see more in Ref. [18]). The loops in all 
three structures have the same dimensions, thus leading to the same magnetic field 
period fioAH = 1.24 mT. The strips forming the structures are 0.13 iim wide and 
the film thickness t — 34 nm. In all the experimental data we show, the parabolic 
background (Eq. (11)) is already subtracted in order to allow for a direct compar- 
ison with the theory. In the temperature interval where we measured the T C (H) 
boundary, the coherence length £ (T) is considerably larger than the width w of the 
strips. This makes it possible to use the one-dimensional models for the calculation 
of E LLL (H) and thus T C (H). The basic idea is to consider \^f s \ = constant across 
the strips forming the network and to allow a variation of \^ s \ only along the strips. 
In the simplest approach \ty s \ is assumed to be spatially constant (London limit 
(LL)) [19,20], in contrast to the de Gennes-Alexander (dGA) approach [12,21,22], 
where \^f s \ is allowed to vary along the strips. In the latter approach one imposes: 



at the points where the current paths join. The summation is taken over all strips 
connected to the junction point. Here, x is the coordinate defining the position 
on the strips, and An is the component of the vector potential along x. Eq. (15) 
is often called the generalized first Kirchhoff law, ensuring current conservation 




(15) 



[22]. The second Kirchhoff law for voltages in normal circuits is now replaced by 
the fluxoid quantization requirement (Eq. (14)), which should be fulfilled for each 
closed contour around a loop. 



In Figs. 10-12 the phase boundaries of the three structures are shown. The 
dashed lines are the phase boundaries obtained with the LL, while the solid lines 
give the results from the dGA approach. As we discussed in Section II. B for a 
mesoscopic loop, attaching contacts modifies the confinement topology, so that the 
amplitude of the local LP oscillations is lowered at low magnetic fields. Here as 
well, the inclusion of the leads reduces the amplitude of the oscillations. The dash- 
dotted line in Figs. 10-12 give the result of the dGA calculation where the presence 
of the leads has been included. The values for £(0) obtained from the fits agree 
within a few percent with the £(0) values found independently from the monotonic 
background of T c ($) (see Eq. (11)). 
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FIGURE 10. Experimental T c ($) data for the bola with the parabolic background (Eq. (11)) 
subtracted. The dots are the experimental data points, while the lines correspond to the different 
theoretical results as explained in the text, a) Single period of T c ($), b) A few periods of the 
experimental T c (<!>) curve. 



1.294 


■ 1 1 — i — 

(a) 


1 ' 1 — 




1.296 








~ 1.298 


/' ' 

In \ 
n 




\ 


\-° 

1.300 








1.302 


I 


II 

I . I 


I > 

i 




0.0 0.2 0.4 0.6 0.8 
®AD Q -n 

FIGURE 11. Experimental T c ($) data for the double loop with the parabolic background 
(Eq. (11)) subtracted. The dots are the experimental data points, while the lines correspond 
to the different theoretical results as explained in the text, a) Single period of T c ($), b) A few 
periods of the experimental T c ($) curve. 



First, in Fig. 10, we consider the mesoscopic "bola" - two loops connected by a 
wire. Fink et al. [22] showed that, in the complete magnetic flux interval, the spa- 
tially symmetric solution, with equal orientation of the supercurrents in both loops, 
has a lower energy than the antisymmetric solution. Coming back to the similarity 
between a mesoscopic loop and a hydrogen atom, we discussed in Section II. B, 
we can then compare the bola with a H 2 molecule, where the symmetric and the 
antisymmetric solutions correspond to singlet and triplet states, respectively. In 
fact, T C (<I>) of the bola is the same as for a single loop provided that the length 
of the strip connecting the two loops is short, as confirmed by the LP oscillations 
observed in the experimental T C (<I>) (Fig. 10). 

In what follows we will focus on the phase boundaries of the double (Fig. 11) 
and the triple loop (Fig. 12). To facilitate the discussion we divide the flux period 
in two intervals: flux regime I for $/$ < g or $/$o > (1 — fi 1 ) and flux regime II 
for g < §/Q < (1 — g). In the flux regime I the phase boundaries, predicted 
by the different models, are nearly identical. Near $/$o = 1/2 (flux regime II), 
however, clear differences are found between the dGA approach and the LL. The 
dGA result fits better the experimental data with respect to the crossover point g 
between regimes I and II, and the amplitude of the T c oscillations. Using the dGA 
approach we have calculated the spatial modulation of |\& s | and the supercurrents 
for different values at the T c (<3>) boundary. In the flux regime I |\& s | varies only 
slightly and therefore the results of the LL and the dGA models nearly coincide. 
The elementary loops have an equal fluxoid quantum number (and consequently an 
equal supercurrent orientation) for both the double and the triple loop geometry. 
For the double loop this leads to a cancellation of the supercurrent in the middle 
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FIGURE 12. Experimental T c ($) data for the triple loop with the parabolic background 
(Eq. (11)) subtracted. The dots are the experimental data points, while the lines correspond 
to the different theoretical results as explained in the text, a) Single period of T c ($), b) A few 
periods of the experimental T c (§) curve. 



strip, while for the triple loop structure the fluxoid quantization condition (Eq. (14)) 
results in a different value for the supercurrent in the inner and the outer loops. 
As a result, the common strips of the triple loop structure carry a finite current. 

In the flux regime II, qualitatively different states are obtained from the LL and 
the dGA approach: the states calculated within the dGA approach have strongly 
modulated \^ s \ along the strips. This is most severe for the double loop: \I/ S 
shows a node (|\I/ S | = 0) in the center of the common strip, the phase if having a 
discontinuity of n at this point. This node is a one-dimensional analog of the core of 
an Abrikosov vortex, where the order parameter also vanishes and the phase shows 
a discontinuity. In Fig. 13 the spatial variation of \^ s \ along the strips is shown 
for $/$ = 0.36 close to the crossover point g. The dashed curve gives \^ s \ in flux 
regime I, which is quasi-constant. The strongly modulated solution, which goes 
through zero in the center, is indicated by the solid line. Although there exists a 
finite phase difference across the junction points of the middle strip, no supercurrent 
can flow through the strip due to the presence of the node. This node is predicted 
to persist when moving below the phase boundary into the superconducting state 
[23,24]. Already in 1964 Parks [25] anticipated that, in a double loop, "a part of the 
middle link will revert to the normal phase", and that "this in effect will convert 
the double loop to a single loop" , giving an intuitive explanation for the maximum 
in T c ($) at $/$ — 1/2- Such a modulation of \^ s \ is obviously excluded in the 
LL, where the loop currents have an opposite orientation and add up in the central 
strip, thus giving rise to a rather high kinetic energy. An extra argument in favor 
of the presence of the node is given by the much better agreement for the crossover 
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FIGURE 13. Calculated variation of \^ s \ along the circumference of the double loop, at the 
phase boundary ($/$o = 0.36) . The dashed line is the solution with | *f? a | nearly spatially constant, 
while the solid line is the state with a node in the center of the strip connecting points A and D. 

point g when the presence of the leads is taken into account in the calculations (see 
dash-dotted line in Fig. 11). 

For the triple loop structure the modulation of |^ s | is still considerable in flux 
regime II, but it does not show any nodes. Therefore the supercurrent orientations 
can be found from the fluxoid quantum numbers {fii}, obtained from integrating 
the phase gradients along each individual loop. When passing through the crossover 
point between flux regime I and regime II only the supercurrent in the middle loop 
is reversed, while increasing the flux above $/$o = 1/2 implies a reversal of the 
supercurrent in all loops. 

Surprisingly, the behavior of a microladder with a linear arrangement of m loops 
appears to be qualitatively different for even and for odd m in the sense that m 
determines the presence or absence of nodes in the common strips. For an infinitely 
long microladder \ty s \ was found to be spatially constant below a certain $ < $ c 
[26], which is analogous to the states we find in flux regime I. For fluxes $ > $ c 
modulated \^! s \ states, with an incommensurate fluxoid pattern, were found. At 
$/<E»o = 1/2, nodes appear at the center of every second common (transverse) 
branch. 

B 2D clusters of antidots 

As a 2D intermediate structure between individual elements A and their huge 
arrays (Fig. 3), we have studied the superconducting micro-square with a 2x2 
antidot cluster. In this case A= "antidot". The micro-square with the 2x2 antidot 




cluster consists of a 2x2 /im 2 superconducting square with four antidots (i.e. square 
holes of 0.53 x 0.53 fim 2 ). A Pb/Cu bilayer with 50 nm of Pb and 17 nm of Cu 
was used as the superconducting film for the fabrication of this structure. The 
thin Cu layer was deposited on the Pb to protect it from oxidation and to provide 
a good contact-layer for the wire-bounding to the experimental apparatus. An 
AFM image of the Pb/Cu 2x2 antidot cluster, is shown in Fig. 14 together with 
a reference sample (i.e. a Pb/Cu micro-square of 2x2 fim 2 without antidots). 
The Pb(50 nm)/Cu(50 nm) bilayer behaves as a Type II superconductor with a 
T c0 = 6.05 K, a coherence length, £(0) ~ 35 nm and a dirty limit penetration 
depth, A(0) 76 nm. The T C (H) measurements on the reference sample revealed 
characteristic features originating from the confinement of the superconducting 
condensate by the dot geometry (see Section II. C). The additional features observed 
in the T C (H) phase boundary of the antidot cluster can therefore be attributed to 
the presence of the antidots. 

The experimental T C (H) phase boundary is shown in Fig. 15. It was measured 
by keeping the sample resistance at 10% of its normal state value and varying the 
magnetic field and temperature [27]. Strong oscillations are observed with a peri- 
odicity of 26 G and in each of these periods, smaller dips appear at approximately 
7.5 G, 13 G and 18 G. The parabolic background superimposed on T C (H) can again 
be described by Eq. (11). 

Defining a flux quantum per antidot as $o = h/2e = BS, where B = fi H and 
S is an effective area per antidot cell (S— 0.8 pm 2 ), the minima observed in the 
magnetoresistance and the T C (H) phase boundary at integer multiples of 26 G can 
be correlated with a magnetic flux quantum per antidot cell, $ = n$ - The ones 
observed at 7.5 G, 13 G and 18 G correspond to the values $/<I>o=0.3, 0.5 and 0.7. 
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FIGURE 15. Experimental phase boundary, AT C (H) for the Pb/Cu 2x2 antidot cluster. 
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FIGURE 16. (a) Theoretical phase boundary, T c (&/$ a ), calculated in the London limit of 
the Ginzburg-Landau theory without any fitting parameter (solid line). All possible parabolic 



The solutions obtained from the London model (LL) define a phase bound- 
ary which is periodic in $ with a periodicity of $o- Within each parabola 
AT C = 7($/$o) 2 > where the coefficient 7 characterizes the effective flux pene- 
tration trough the unit cell. The 7-value is determined by the combination of A 
and the effective size of the current loops. In Fig. 16, the first period of this phase 
boundary, AT C (<3>) = T c0 - T c ($) versus $/$ is shown. Note that there are six 
parabolic solutions given by a different set of flux quantum numbers {rii}, each 
one defining a specific vortex configuration. In Fig. 16a, this is indicated by the 
numbers shown inside the schematic drawings of the antidot cluster. Note that 
some vortex configurations are degenerate. 

From all these possible solutions, for each particular value of $/$o, only the 
branch with a minimum value of AT C ($) is stable (indicated with a solid line in 
Fig. 16a). Note that for the phase boundary, calculated within the ID model of 4 
equivalent and properly attached squares, no fitting parameters were used since the 
variation of T C (<I>) was calculated from the known values for £ and the size. One 
period of the phase boundary of the antidot cluster is composed of five branches 
and in each branch a different stable vortex configuration is permitted. For the 
middle branch (0.37< $/$o <0.63), the stable configuration is the diagonal vortex 
configuration (antidots with equal rii at the diagonals) instead of the parallel state 
(dashed line in Fig. 16a). 

The net supercurrent density distribution circulating in the antidot cluster for 
different values of $/<E»o has been determined using the same approach. Circular 
currents flow around each antidot. For the states rij=0 and rii—1 currents flow 
in the opposite direction, since currents corresponding to nj=0 must screen the 
flux to fulfil the fluxoid quantization condition (Eq. (14)), whereas for nj=l they 
have to generate flux. At low values of $/$o> currents are canceled in the internal 
strips and screening currents only flow around the cluster. When we enter the field 
range corresponding to the second branch of the phase boundary, a vortex (rij=l) 
is pinned around one antidot of the cluster (see Fig. 16a). At the third branch, the 
second vortex enters the structure and is localized in the diagonal. In the fourth 
branch of the phase boundary the third vortex is pinned in the antidot cluster. 
And finally, the current distribution for the fifth branch is similar to that of the 
first branch although currents flow in opposite direction. 

Figure 16c shows the first period of the measured phase boundary T c ($) after 
subtraction of the parabolic background. For all measured samples, the first period 
of the experimental phase boundary is composed of five parabolic branches with 
minima at $/$ = 0, 0.3, 0.5, 0.7, I. If we compare it with the theoretical prediction 
given in Figure 16a, the overall shape can be reproduced although the experimental 
plot has two major peaks at $/$o — 0.2 and 0.8 whereas the theoretical curve only 
predicts cusps around these positions. 

The agreement between the measured and the calculated T c ($) is improved if we 
assume that the coefficient 7 can be considered as a fitting parameter. This seems 
to be feasible if we take into account the simplicity and limitation of the used ID 
model. Due to the relatively large width of the strips forming the 2x2 cluster, 



the sizes of the current loops can change since they are " soft" in this case and not 
defined very precisely. 

As a result, the coefficient 7 can not be treated as a known constant. If we use 
it as a free parameter (Fig. 16b) then the curvature of all parabolae forming T C (H) 
can be changed and the calculated T C (H) curve becomes closer to the experimental 
one though the amplitude of the maxima at $/$o = 0.2 and 0.8 is still lower than in 
the experiment (Fig. 16c). The discrepancy in the description of the amplitude of 
the maxima at $/$o= 0.2 and 0.8 could also be related to the pinning of vortices by 
the antidot cluster when potential barriers between different vortex configurations 
may appear. At the same time, the achieved agreement between the positions 
of the measured and calculated minima of the T C (H) curves confirms that the 
observed effects are due to fluxoid quantization and formation of certain stable 
vortex configurations at the antidots. 



IV SUPERCONDUCTING FILMS WITH AN 

ANTIDOT LATTICE 



Laterally nanostructured superconducting films having regular arrays of antidots 
are convenient model objects to study the effects of the confinement topology on the 
T C (H) phase boundary in two different regimes. The first (or "collective") regime 
corresponds to the situation where all elements A, forming an array, are coupled. 

From the experimental T C (H) data on antidot clusters we expect for films with 
an antidot lattice, higher critical fields at $ = n<3> , which is in agreement with 
the appearance of the T C (H) cusps at $ = n$ m superconducting networks [28]. 
Here, the flux $ is calculated per unit cell of the antidot lattice. 

On the other hand, by applying sufficiently high magnetic fields, the individual 
circular currents flowing around antidots, can be decoupled and the crossover to 
a "single object" behaviour could be observed. In this case the relevant area for 
the flux is the area of the antidot itself and we deal with surface superconductivity 
around an antidot. 

Figure 17 shows the critical field for a Pb(50 nm) sample with a square antidot 
lattice (period d — 1 fiin and the antidot size a = 0.4 jum). The T C (H) boundary is 
determined at 10% of the normal state resistance, R n . In this graph two distinct 
periodicities are present. 

Below ~ 8 mT cusps are found with a period of 2.07 mT, corresponding to one 
flux-quantum per lattice cell. These cusps or "collective oscillations" [29] are remi- 
niscent of superconducting wire networks [28] and arise from the phase correlations 
between the different loops which constitute the network. 
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FIGURE 17. Critical field of a Pb(50nm) film measured at 10% R n , with d — lfrni, a = 0.4/^m. 
The inset shows a zoom of the first cusp at different criteria 50% R n , 10% R n and 0.5% R n . 



These cusps are obtained by narrowing the minima at n$o with increasing size N 
of the N x N antidot cluster (see the sharpening of the minima at $ = 0, $o m 
Fig. 18; note that the phase boundary in the N — > oo case has a similar shape as the 
lowest energy level of the Hofstadter butterfly [2,31]). An important observation 
is that the amplitude of these "collective" oscillations depends upon the choice of 
the resistive criterion. This is similar to the case of Josephson junction arrays and 
weakly coupled wire networks [30] where phase fluctuations dominate the resistive 
behavior. The inset of Fig. 17 shows the first collective period, measured using 
three different criteria. As the criterion is lowered the cusps become sharper and 
the amplitude increases well above the prediction based on the mean field theory 
for strongly coupled wire networks [28]. At the same time, cusps appear at rational 
fields $/<l > o=l/4, 1/3, 1/2, 2/3 and 3/4 arising from the commensurability of the 
vortex structure with the underlying lattice. 

Above ~ 8 mT, the collective oscillations die out and "single object" cusps ap- 
pear, having a periodicity which roughly corresponds to one flux quantum $ per 
antidot area, a 2 . These cusps are due to the transition between localized super- 
conducting edge states [29] having a different angular momentum L. These states 
are formed around the antidots and are described by the same orbital momentum 
introduced in Section II. C for dots. 
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FIGURE 18. Calculations of the first T C (H) period for an N x N antidot system 
(N = 1, 2, 3, 4, oo ) in the London limit. 
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FIGURE 19. The critical field of Fig. 17 normalized by <J> /27r£ 2 (T) versus the applied field. 
The dashed line (right axis) shows the theoretical result for a single circular hole with a radius 
R = 0.24 nm. 

Figure 19 shows the same critical field as shown in Fig. 17, but normalized by 



the upper critical field H c2 of a plain film without antidots, HqH c2 = $ /2vr^ 2 (T) 
( £(0) = 36 nm). The dashed line is the calculation of the reduced critical field 
for a plain film with a single circular antidot with radius R = 0.24 jum. The 
positions of the cusps correspond reasonably well to the experimental ones, taking 
into account that the model only considers a single hole. From this comparison, 
an effective area irR 2 = 0.187 /zm 2 is determined which is close to the experimental 
value a 2 = 0.16 /im 2 . 

From Figures 17 and 19 it is possible to show that the transition from the network 
regime to the "single object" regime takes place at a temperature T* given by the 
relation w ~ 1.6 £(T*), (where w is the width of the superconducting region between 
two adjacent antidots). 

Experiments on systems with other antidot sizes demonstrate that the a/d ra- 
tio determines the relative importance of the "collective regime" and changes the 
cross-over temperature T*. The relation w ~ 1.6 £(T*), seems nevertheless to hold 
reasonably well and is similar to the transition from bulk nucleation of superconduc- 
tivity to surface nucleation in a thin superconducting slab parallel to the magnetic 
field [15], which happens at a temperature T cr satisfying w — 1.8 £(T cr ). 

Comparing the bulk H c2 (T) curve with the T C (H) boundary for films with an 
antidot lattice, we clearly see a qualitative difference between the two, caused by 
the lateral micro-structuring. In the network limit, T C (H) can be related to the 
lowest Em level in the Hofstadter butterfly [2,31] with pronounced cusps at n$ 
and the substructure within each period. Reducing the size of the antidot, we are 
modifying T C (H) substantially but the cusps at n$ are still clearly seen [32]. 

Finally, by introducing antidot lattices we stabilize the novel flux states, which 
can be briefly characterized as follows. For relatively large antidots sharp cusp-like 
magnetization anomalies appear at the matching fields H m . These anomalies are 
caused by the formation of the multi-quanta vortex lattices at each subsequent H m 
[33-36]. The multi-quanta vortex lattices make possible a peaceful coexistence of 
the flux penetration at the antidots and the presence of a substantial superfluid 
density in the space between them. This leads to a very strong enhancement of the 
critical current density in films with an antidot lattice. For smaller antidots the 
vortices are forced to occupy the interstitial positions after the saturation of the 
pinning sites at antidots [36-38] . This leads to the formation of the novel composite 
flux-line lattices consisting from the interpenetrating sublattices of weakly pinned 
interstitial single-quantum vortices and multi-quanta vortices strongly pinned by 
the antidots. When the interstitial flux-line lattice melts, it forms the interstitial 
flux-liquid coexisting with the flux solid at the antidots. 

V CONCLUSIONS 

We have carried out a systematic study of and quantization phenomena in nanos- 
tructured superconductors. The main idea of this study was to vary the boundary 
conditions for confining the superconducting condensate by taking samples of dif- 



ferent topology and, through that, to modify the lowest Landau level E LLL (H) and 
therefore the critical temperature T C (H). Three different types of samples were 
used: (i) individual nanostructures (lines, loops, dots), (ii) clusters of nanoscopic 
elements - ID clusters of loops and 2D clusters of antidots, and (iii) films with 
huge regular arrays of antidots (antidot lattices). We have shown that in all these 
structures, the phase boundary T C (H) changes dramatically when the confinement 
topology for the superconducting condensate is varied. The induced T C (H) vari- 
ation is very well decribed by the calculations of Elll(H) taking into account 
the imposed boundary conditions. These results convincingly demonstrate that 
the phase boundaryin T C (H) of nanostructured superconductors differs drastically 
from that of corresponding bulk materials. Moreover, since, for a known geometry 
Elll{H) can be calculated a priori, the superconducting critical parameters, i.e. 
T C (H), can be controlled by designing a proper confinement geometry. While the 
optimization of the superconducting critical parameters has been done mostly by 
looking for different materials, we now have a unique alternative - to improve the 
superconducting critical parameters of the same material through the optimiza- 
tion of the confinement topology for the superconducting condensate and for the 
penetrating magnetic flux. 
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